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ABSTRACT 

We present black hole mass-bulge velocity dispersion relation for a complete sample of 75 
soft X-ray selected AGNs: 43 broad line Seyfert Is and 32 narrow line Seyfert Is. We use 
luminosity and FWHM(H/3) as surrogates for black hole mass and FWHMQOIII]) as a surrogate 
for the bulge velocity dispersion. We find that NLSls lie below the M BK -a relation of BLSls, 
confirming the Mathur et al. (2001) result. The statistical result is robust and not due to any 
systematic measurement error. This has important consequences towards our understanding of 
black hole formation and growth: black holes grow by accretion in well formed bulges, possibly 
after a major merger. As they grow, they get closer to the Mbh _ c relation for normal galaxies. 
The accretion is highest in the beginning and dwindles as the time goes by. Our result does not 
support theories of Mbh-c relation in which the black hole mass is a constant fraction of the 
bulge mass/ velocity dispersion at all times or those in which bulge growth is controlled by AGN 
feedback. 
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Active galaxies are "active" because they ac- 
crete matter on to the supermassive black holes. 
However, whether the accretion leads to significant 
growth of the nuclear black hole has been a matter 
of some debate. New results on X-ray background 
and the better determination of local black hole 
mass density have led to the conclusion that in- 
deed, black holes grow during their active phase 
(e.g. Barger et al. (2001); Aller & Richstone 
(2002); Yu & Tremainc (2002)). 

The above result needs to be understood in the 
context of the observed tight correlation between 
the central black hole Mbh and the stellar ve- 
locity dispersion er* of the bulge in a galaxy (e.g. 
Gebhardt et al. (2000a); Ferrarese & Merritt 
(2000); Merritt & Ferrarese (2001)). Measuring 
black hole masses using stellar and gas dynamics 
(Kormendy & Richstone (1995) and references 
therein), these authors found that log Mbh = 



a + frxlog (er*/ero) with A/bh in units of M q and 
(To =200 km s . The slopes b of this relation vary 
between b — 3.75 (Gebhardt et al. 2000a) and 
b = 5.27 (Ferrarese & Merritt 2000). Most re- 
cently, Tremaine et al. (2003) estimated b = 4.02 
and a — 8.13 for a sample of 31 nearby galaxies. 
Interestingly, the above relation for normal galax- 
ies also extends to active galaxies (Gebhardt et al. 
2000b; Ferrarese et al. 2001). 

The above two results imply that the formation 
and growth of the nuclear black hole and the bulge 
in a galaxy are intimately related, and several the- 
oretical models have attempted to explain the ob- 
served Mbh - 0- relation (e.g. Haenhnelt (2003); 
Haehnelt et al. (1998); Adams et al. (2001) and 
King (2003)). To understand the link between 
the black hole and the bulge, it is important to 
determine whether (a) black hole mass is a con- 
stant fraction of the bulge mass, or bulge velocity 
dispersion, at all times, or (b) during some accret- 
ing phase, the Mbh - o~ relation is not followed by 
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AGNs. The former can be obtained if, e.g., the 
growth of a black hole matches the growth of its 
surrounding bulge exactly during a merger. Alter- 
natively, feedback from a black hole can limit the 
bulge growth. In case (b), the accretion history of 
a black hole is not tied to the bulge growth. It is 
still possible that today's highly accreting AGNs 
do lie on the Mbh~ & relation for normal galaxies, 
and may eventually leave the relation as their BHs 
grow. Or, today's highly accreting AGNs may lie 
below the Mbh~ cr relation for normal galaxies, 
and would eventually reach the relation as the ac- 
tive BHs become dead. Indeed, in a recent model 
by Miralda-Escude & Kollmeier (2003), which ex- 
plicitly couples accretion with the stellar system 
in the bulge, the observed Mbh _ cr relation is the 
final relation at the end of the accretion process. 

It is of interest, therefore, to follow the tracks 
of AGNs on the Mbh - cr plane. Since black holes 
would grow fastest with high accretion rates, ac- 
tive galaxies with close to Eddington accretion are 
perhaps the best candidates. At low rcdshift, a 
lot of observational evidence suggests that nar- 
row line Seyfert 1 galaxies (NLSls; a subclass of 
Seyfert galaxies with full width at half maximum 
of H/3 lines less than 2000 km s _1 (Osterbrock & 
Pogge 1985)) accrete at close to Eddington rate 
(e.g. Pounds et al. (1995); Grupe (2004) and 
references there in). Mathur ct al. (2001) ar- 
gued that NLSls do not follow the Mbh - cr rela- 
tion. This, however, was only a suggestive result 
because the sample size was small and neither the 
black hole masses nor the velocity dispersions were 
accurately measured. Wandel (2002) confirmed 
the Mathur et al. result with a sample of 55 AGNs 
(see also Bian & Zhao (2003)). On the other hand, 
Wang & Lu (2001) have argued that both BLSls 
and NLSls follow the Mbh - cr relation. Clearly, 
it is important to test whether NLSls occupy a 
distinct region in the Mbh _ cr plane using a large, 
homogenous sample. In this paper we present our 
results based on a complete sample of 110 soft X- 
ray selected AGN. 

Note also that NLSls are interesting objects as 
they occupy one extreme end of the "eigenvector 
1" relation of AGNs (Boroson & Green (1992)). 
The most widely accepted paradigm for NLSls is 
that they accrete at close to Eddington rate and 
have smaller black hole masses for a given lumi- 
nosity compared to BLSls. Finding their locus on 



the Mbh _ ct plane is therefore a worthwhile exper- 
iment anyway as we will either find that occupy 
a distinct region compared to BLSls or that they 
don't. The first option is interesting for the rea- 
sons discussed above. On the other hand if we 
find that NLSls follow the Mbh _ ct relation like 
the BLSls, it has important implications towards 
our understanding of the AGN phenomenon. As 
noted above, we already have a good evidence for 
smaller BH masses of NLSls, at a fixed luminosity. 
If they follow the Mbh _ ct relation, it would imply 
that NLSls preferentially reside in galaxies with 
bulges of smaller velocity dispersion. This would 
be a direct evidence for dependence of AGN prop- 
erties on their large scale galactic environment. 

In section 2 we discuss the sample selection 
and the methodology to determine the black hole 
masses using the widths of the H/3 lines. Widths of 
the forbidden [OIII] lines are used as surrogates for 
bulge velocity dispersion. In section 3 we present 
the results and the discussion is in section 4. 

2. Estimates of Mbh and a 

2.1. The Sample 

The sample contains all bright soft X-ray se- 
lected AGN from the ROSAT All-Sky Survey 
(RASS, Voges et al. (1999)). The selection crite- 
ria are given in Thomas et al. (1998) and Grupe 
ct al. (1999). About half of these sources are 
NLSls (51 objects) and 59 are broad line seyfert 
Is (BLSls). NLSls and BLSls show similar distri- 
bution in their rcdshifts, luminosities and equiv- 
alent widths of H/3 (Grupe et al. 2004). From 
this original sample of 110 AGNs, we removed all 
the objects in which [OIII] lines were severely un- 
resolved leading to errors in the FWHM measure- 
ments, with S/N< 3. This left us with a sample 
of 75 AGNs, 32 NLSls and 43 BLSls. 

2.2. Black hole mass 

Based on the reverberation mapping of the H/3 
line of 28 PG quasars, Kaspi et al. (2000) 
derived an empirical relation between the width 
of the H/3 line and the central black hole mass. 
They found that log M B h =5.17+log i? B LR + 
2 x (logFWHM(H/3) - 3) with M BH given in 
units of M Q and FWHM (H/3) given in units of 
km s _1 . i?BLR is the radius of the broad emission 
line region (BLR) and is larger for more luminous 
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sources: log .Rblr = 1.52 + 0.70 x (logAi 5 i 00 - 37) 
where -Rblr given in units of light days, L5100 
is the monochromatic luminosity at 5 100 A, and 
Aisloo is iu units of Watts. This relation is well 
calibrated, albeit with some scatter, and can be 
used to estimate Mbh of active galaxies using 
FWHM(H/3) and luminosity. Many authors have 
used it to derive Mbh of Seyfert galaxies and 
quasars (e.g. Laor (1998) and Wandel (1999)) 
and we do the same (except for NGC 4051, which 
is known not to follow the radius luminosity re- 
lation, and so we use the BH mass as measured 
by reverberation (Peterson et al. 2000)). For our 
sample, FWHM(H/3) and optical luminosities are 
given in Grupe et al. (2004). 

2.3. Velocity dispersion 

The stellar velocity dispersion cr* in the bulge 
of a galaxy can be measured by the widths of the 
Call triplet absorption features at 8498.0, 8542.1, 
and 8662.1 A. In a sample of 85 AGN, Nelson & 
Whittle (1995, 1996) found a moderately strong 
correlation between cr* and FWHM([OIII]), the 
full width at half maximum of the [OIII]A5007 
line. Therefore, Nelson (2000) suggested that 
the FWHMQOIII]) can be used as a surrogate of 
cr* with cr* w £7 [0 iii] = FWHM([OIII])/2.35. 
This result was confirmed by Boroson (2003), 
who stated that the [OIII] width can predict the 
black hole mass to a factor of 5 assuming the 
Mbh~ cr relation. Shields et al. (2003) used the 
FWHM[(OIII)] as a surrogate of the bulge stellar 
velocity dispersion cr* for a large sample of AGN 
and concluded that it can be extended even to red- 
shifts as high as z«3. We also use FWHMQOIII]), 
given in Grupe et al. (2004), as an estimate of cr* 
for our sample. 

3. Results 

Fig. 1 plots the black hole mass derived from 
FWHM(H/3) vs. velocity dispersion cr derived 
from the FWHMQOIII]) for our sample of 75 
AGNs. The solid line is the relation of Trcmaine et 
al. (2003) obtained for normal galaxies. Clearly 
there is a large scatter around this relation, and 
most likely implies that the surrogates do not re- 
produce Mbh and cr accurately. A comparison 
of the black hole masses derived from our data 
with reverberation mapping results (Kaspi et al. 



(2000); Peterson et al. (2000); Wandel (1999); 
Kollatschny (2003)) shows that the black hole 
masses agree on average, with a random scatter 
of about 0.2 dex. Considering the scatter in the 
radius -luminosity relation, variability and the un- 
known geometry of the broad line region, we con- 
servatively estimate the error on BH mass mea- 
surement to be 0.5 dex. 

The observational error in cr depends on the 
strength of the [OIII] line. As shown in Grupe 
(2004) there is a strong anti-correlation between 
the FWHMQOIII]) and equivalent width of [OIII]. 
Plus the weaker the [OIII] line the stronger the 
Fell emission becomes (e.g. Boroson & Green 
(1992); Grupe (2004)) which makes the line mea- 
surements in the objects with broader [OIII] more 
uncertain than those with narrow [OIII] emis- 
sion. The errors of the FWHMQOIII]) are given 
in Grupe et al. (2004) and are of the order of 0.2 
dex. Clearly, errors on both quantities, Mbh and 
a are large for individual objects. However, we are 
interested in statistical differences in the two pop- 
ulations, of a large number of BLSls and NLSls, 
and not on exact values of individual sources. 

Figure 1 shows that BLSls and NLSls occupy 
two distinct regions in the Mbh~ cr plane. For 
a given velocity dispersion NLSls tend to show 
smaller smaller black hole masses than BLSls. If 
true, this is an important result. However, be- 
fore coming to that conclusion, we have to make 
sure that the result is not spurious. A spurious 
result may be obtained if the black hole masses 
of NLSls are systematically underestimated or if 
their velocity dispersions are systematically over- 
estimated relative to BLSls. We will check these 
two cases below. 

Are the black hole masses of NLSls wrong? 
Fig. 2 shows the cumulative fraction of the dis- 
tributions of the inferred black hole masses of 
NLS1 (solid line) and BLSls (dashed line) for a 
Kolmogorov-Smirnov (KS) test. The plot clearly 
shows that NLS1 and BLSls have different distri- 
butions of the black hole masses. In general, more 
luminous AGNs have higher black hole masses, 
but for a given luminosity NLSls have black hole 
masses about an order of magnitude lower than 
BLSls (note that the BLSls and NLSls in our 
sample have similar luminosities; §2.1). This re- 
sult confirms earlier findings of, e.g. Wandel & 
Boiler (1998) and Peterson et al. (2000) and is 



3 



unlikely to be spurious. In fact, NLSls have nar- 
rower broad emission lines because of the smaller 
black hole masses. If they had masses similar to 
the BLSls, their BLRs would have to be relatively 
farther away from the black hole. This, however, 
is not the case; Peterson et al. (2000) have found 
that NLSls and BLSls follow the same relation 
between the BLR size and luminosity. We thus 
conclude that there is a real difference in the black 
hole mass distribution of BLSls and NLSls. 

Are the estimates of velocity dispersion wrong? 
If FWHMQOIII]) is not a good surrogate, our 
estimates of a may well be wrong. This could 
produce dichotomy between NLSls and BLSls if 
a [Oiii] ~ °~* is systematically different for the two 
classes. The first indication that this is not the 
case comes from the similarity of distribution of 
their a (Fig. 3). A KS test shows that the two 
classes do not show any significant difference. 

One might overestimate the [OIII] widths if the 
spectral resolution is low and lines are not re- 
solved. We correct for the instrumental line broad- 
ening in our measurements (Grupe et al. 2004). 
Moreover, if resolution was a problem, we would 
have seen a clustering of line widths close to the in- 
stumcntal line widths. Instead, we observe a wide 
range of widths, and the distribution of widths for 
BLSls and NLSls is similar. This implies that any 
problem with resolution is not artificially increas- 
ing the line widths of NLSls only. 

One problem which might affect NLSls only is 
the strong Fell emission close to O[III]A5007 emis- 
sion line. For our entire sample, Fell contribution 
has been subracted before making measurements 
on the [OIII] lines (Grupe et al. 2004). If Fell 
contribution was systematically undersubtracted, 
it will lead to overestimation of [OIII] line widths. 
With this in mind, we re-examined the Fell sub- 
tracted spectra of all the NLSls in our sample. 
We found that Fell emission was not undersub- 
tracted, and may even be slightly oversubtracted, 
in all but but two cases. Deleting these two ob- 
jects with poorer S/N from our sample does not 
change the statistical result. 

A more detailed look at many [OIII] lines shows 
that they have blue asymmetry in a number of 
AGNs (e.g. Goncalves et al. (1999); Leighly 
(2000)), possibly resulting from outflows. These 
observations need further scrutiny since the asym- 
metry might have resulted in erronously large 



measurements of FWHM([OIII]). (Note, however, 
that a part of the blue asymmetry may be a result 
of oversubtracting the Fell emission from the red 
part of the line.) To correct for this possible prob- 
lem, we re-measured the width of the [OIII] line 
as 2 x half width at half maximum of the red part 
of emission line. While the new measurements re- 
duce individual values of <r, the BLSls and NLSls 
still ccupy distinct regions in the Mbh~ o plot. 
This is clearly seen from figure 4 which plots the 
cumulative distribution of the ratio Mbh/c for 
BLSls and NLSls. The two classes are signifi- 
cantly different, with formal K-S test probability 
of being drawn from the same population < 0.001. 

4. Discussion 

In this paper we present distribution of black 
hole masses and bulge velocity dispersions for a 
complc sample of AGNs. Black hole masses are 
estimated from a well calibrated relation between 
the width of the H/3 emission line, luminosity and 
A^bh- Bulge velocity dispersion is estimated from 
the width of the narrow [OIII] emission line. Nei- 
ther of these are accuarate measurements of the 
two quantities and errors on the values of Mbh 
and a for individual objects are large. The results 
presented here are therefore statistical, and com- 
pare the broad distributions of the two classes of 
AGNs, the BLSls and the NLSls. 

We find that BLSls and NLSls occupy two dis- 
tinct regions in the Mbh - cr plane. This does not 
appear to be a result of systematically underesti- 
mating black hole masses or systematically over- 
estimating the [OIII] line widths of NLSls. This 
might be still a spurious result if the narrow emis- 
sion line region of NLSls is somehow much closer 
to the center of the galaxy and so does not trace 
the bulge velocity dispersion. While we cannot 
rule out this possibility, there are no observations 
supporting this fact either. 

After carefully investigating all the options, we 
come to the conclusion that the result is unlikely 
to be spurious in that NLSls and BLSls do have 
different ratios of black hole mass to bulge ve- 
locity dispersion. So if BLSls follow the Mbr— 
cr relatation for normal galaxies, then NLSls do 
not. Needless to say, this result will have to be 
confirmed with accurate measurements of black 
hole masses and bulge velocity dispertions. If con- 
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firmed, it has important consequences towards our 
understanding of black hole formation and growth. 
We find that the black hole mass is not a constant 
fraction of the bulge mass or of bulge velocity dis- 
persion at all times. In other words, growth of a 
black hole by accretion does not match the growth 
of its surrounding bulge exactly or the accretion 
process itself does not increase the bulge veloc- 
ity dispersion. Thus, our results are also inconsis- 
tent with the models of M^n-a relation in which 
feedback from a black hole controls the growth 
of the bulge. Our resuls support a scenario in 
which black holes grow by accretion in well formed 
bulges, possibly after a major merger. Perhaps, 
the surge in accretion is a result of that merger 
itself. As they grow, they get closer to the Mbh - 
a relation for normal galaxies. The accretion rate 
is highest in the beginning and dwindles as the 
time goes by. This scenario is consistent with 
the recent theoretical model of Miralda-Escude & 
Kollmcicr (2003) in which the Mbh-c relation 
is the final state at the end of the accretion pro- 
cess. The observation that the broad line Seyfert 
Is lie close to the Mbh-c relation then tells us 
that the black hole mass growth at low accretion 
rate is not significant. This scenario is also consis- 
tent with the proposal of Grupe et al. (1999) and 
Mathur (2000) that NLSls are younger members 
of the Seyfert population. Presumably, the insight 
we have got from studying the Seyfert galaxies is 
also applicable to quasars. To understand further 
the history of BH and bulge growth, similar ob- 
servations at high redshift would be valuable. 

We thank Jordi Miralda-Escude, Brad Peter- 
son, Rick Pogge, Greg Shields, David Weinberg 
and Bev Wills for helpful discussions. This re- 
search has made use of the NASA/IPAC Extra- 
galactic Database (NED) which is operated by the 
Jet Propulsion Laboratory, Caltcch, under con- 
tract with the National Aeronautics and Space 
Administration. 
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Fig. 1. — Velocity dispersion <7[om] vs. log 
Mbh(H/3). NLSls are marked as open circles and 
BLSls as filled squares. Black hole masses are 
given in units of Mq. The solid line marks the 
relation of Tremaine et al. (2003). The cross at 
the bottom right hand corner represents a typical 
error bar. 
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Fig. 3. — Cumulative fraction of a KS test of the 
distributions of the stellar velocity dispersion a 
given in units of km s _1 . The distribution of 
NLSls is shown as a solid line and BLSls are 
shown as a dashed line. 




Fig. 2. — Cumulative fraction of a KS test of the 
black hole mass distributions of NLSls (solid line) 
and BLSls (dashed line) given in units of M Q . 
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Fig. 4. — Cumulative fraction of a KS test of the 
distributions of the black hole mass Mbh divided 
by the stellar velocity dispersion a. The distribu- 
tion of NLSls is shown as a solid line and BLSls 
are shown as a dashed line. 
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